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I "Abstract 

{S) One significant difference between the atmospheres of stars and exoplanets is the presence of condensed particles (clouds 
kjp^or hazes) in the atmosphere of the latter. In current ID models clouds and hazes are treated in an approximate way 
.by raising the surface albedo, or adopting measured Earth cloud properties. The former method introduces errors to 
<^ jthe modeled spectra of the exoplanet, as clouds shield the lower atmosphere and thus modify the spectral features. The 
.latter method works only for an exact Earth-analog, but it is challenging to extend to other planets. 

The main goal of this paper is to develop a self-consistent microphysical cloud model for ID atmospheric codes, which 
.can reproduce some observed properties of Earth, such as the average albedo, surface temperature, and global energy 
i- r budget. The cloud model is designed to be computationally efficient, simple to implement, and applicable for a wide 
Ph .range of atmospheric parameters for planets in the habitable zone. 

I^j ' We use a ID, cloud-free, radiative-convective, and photochemical equilibrium code originally developed by Kasting, 
Pavlov, Segura, and collaborators as basis for our cloudy atmosphere model. The cloud model is based on models used 
Ojby the meteorology community for Earth's clouds. The free parameters of the model are the relative humidity and 
q .number density of condensation nuclei, and the precipitation efficiency. In a ID model, the cloud coverage cannot be 
p-j 'self-consistently determined, thus we treat it as a free parameter. 

We apply this model to Earth (aerosol number density 100 cm -3 , relative humidity 77 %, liquid cloud fraction 40%, 
, & , 'and ice cloud fraction 25%) and find that a precipitation efficiency of 0.8 is needed to reproduce the albedo, average 
surface temperature and global energy budget of Earth. We perform simulations to determine how the albedo and the 
t-H 'climate of a planet is influenced by the free parameters of the cloud model. We find that the planetary climate is most 
^ .sensitive to changes in the liquid water cloud fraction and precipitation efficiency. 

The advantage of our cloud model is that the cloud height and the droplet sizes are self-consistently calculated, both 
.of which influence the climate and albedo of exoplanets. 
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o: 

"1. Introduction expect that exoplanets in the Habitable Zone (HZ) also 

. . . harbor clouds in their atmospheres. 

Clouds on exoplanets influence the temperature-pressure xhe climatic effects of clouds can be mimicked by in- 
3 ' Pr ° file ^ T_P pr ° file ) ° f the atm0S P here owin S t0 their hi S h creasing the surface albedo (see Kasting et al., 1984, and 
"£j -opacity and therefore influence the borders of the Habit- subsequent work). Although this method can capture the 
03 able Zone (see e.g., Forget and Pierrehumbert, 1997; Selsis climatic effects of clouds and reproduce the observed pressure- 
et al., 2007; Kitzmann et al., 2010). They also influence the temperature structure of Earth, the measured global mean 
detectable spectra by reducing the observable atmospheric annual energy budget of Earth cannot be reproduced (Gold- 
features (see e.g., Kaltenegger et al., 2007; Kitzmann et al., blatt and Zahn le, 2011). Furthermore, the effects of clouds 
2011). Clouds are also a source of uncertainty in problems on the observed spectra is neglected by this approach. To 
such as the Faint Young Sun Paradox (see e.g., Sagan and address this issue, the measured cloud properties of Earth 
Mullen, 1972; Goldblatt and Zahnle, 2011). Earth, Venus, can be used m mo dels assuming that the cloud properties 
and to a small extent, Mars all have clouds, therefore we do not change for different exoplanets (Kitzmann et al., 

2010). An alternative approach is to take the cloud heights 
and droplet sizes from observations, and fit the cloud frac- 
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in agreement with observations (Rossow et al., 2005; War- 
ren and Eastman, 2007). A parameterized cloud model 
is difficult to adopt for exoplanets which are not the ex- 
act duplicate of Earth because the number of cloud layers, 
and other cloud properties (such as the condensing mate- 
rial, height of the layers, droplet sizes, cloud fraction) are 
all unconstrained parameters in the general case. 

Due to the large amount of observational and labo- 
ratory data available for Earth, detailed theoretical cloud 
models have been developed by the meteorological commu- 
nity and are available for ID vertical (see e.g., Srivastava, 
1967), 2D (see e.g., Fan et al., 2007; Shima et al., 2009), 
and 3D general circulation models (see e.g., Stephens, 2005) 
We do not use these models because they are either too 
Earth-centric and study one specific cloud type on a lo- 
cal scale (e.g., cumulus or cirrus clouds only), and/or they 
are too computationally expensive and thus not suitable 
for parameter space exploration. Instead we develop a new 
ID model using basic and well-understood physics also em- 
ployed in these models. 

We develop a new ID microphysical cloud model which 
is suitable for a wide range of planetary and atmospheric 
parameters for exoplanets in the HZ. Our cloud model has 
five free parameters, and does not add a significant com- 
putational overhead to existing ID atmosphere models. 

The structure of the paper is as follows. We describe 
the numerical methods and codes used for the atmosphere 
modeling in Sect. 2. The optical properties of liquid and 
solid water droplets and our microphysical cloud model 
is discussed in Sect. 3. We give the initial conditions of 
the simulations and perform convergence tests in Sect. 4. 
We perform four sets of simulations to study the climatic 
effects of the five free parameters of the cloud model in 
Sect. 5. Finally, we discuss and summarize our results in 
Sects. 6 and 7. 

2. Numerical methods 

2.1. The ID radiative- convective and photochemical code 
We use a radiative-convective climate model (Kasting 
et al., 1993; Pavlov et al., 2000) that is coupled to a photo- 
chemical model (Pavlov and Kasting, 2002; Segura et al., 
2003, 2005, 2010; Kaltenegger et al., 2011). The radiative- 
convective model determines the vertical temperature and 
pressure profiles of the atmosphere in radiative equilibrium 
for a fixed chemical composition. The photochemistry 
code determines the vertical mixing ratio profiles of dif- 
ferent chemical species in chemical equilibrium for a fixed 
temperature and pressure profile. The two modules are 
used iteratively until both radiative and chemical equilib- 
ria are simultaneously reached. 

The photochemistry code was originally developed by 
Kasting et al. (1985). It solves for 55 chemical species 
that are linked by 220 reactions and it uses a reverse-Euler 
method to solve the chemical network (Segura et al., 2010). 

The climate model consists of two parts. The first part 
calculates the transfer of stellar radiation in the visible 



and near-infrared. The second part calculates the trans- 
fer of thermal infrared radiation of the planet. Knowing 
the difference in the incoming and outgoing fluxes and the 
heat capacity of the atmospheric gases, we calculate the 
temperature change of each layer: 

dT = g_dF_ 

dt c p dP' { ' 

where T is the temperature of the layer, dt is the time 
step, g is the gravitational constant, c p is the total heat 
capacity of the layer that depends on the composition of 
the atmosphere, dF is the net flux propagating through 
layer (difference between incoming and outgoing fluxes), 
dP is the pressure difference at the top and bottom of the 
layer. If the temperature gradient (dT/dz - lapse rate) 
in a given layer is larger than the adiabatic lapse rate, 
convective adjustment is performed. 

2.2. Changes in the ID code 

To implement clouds into the ID model, we need high 
vertical resolution (small dz), because clouds only occupy 
a small vertical portion of the atmosphere. Typically the 
height of the atmospheres we consider are about 70 km 
with a cloud layer of a few hundred meters, although deep 
convective clouds around the equator can be several kilo- 
meters thick (Mace et al., 2009). Therefore the number of 
vertical layers is a free parameter in our code and conver- 
gence tests are used to determine the appropriate number 
of layers to be used (see Sect. 4.2). 

We decrease the time step of the calculations and in- 
crease the number of iterations compared to the cloud-free 
version code due to the high resolution of our grid (see 
Sect. 4.1 for more details). We update the convergence 
test of the original code to examine the relative change 
of temperature and mixing ratio of important greenhouse 
gases between two successive calls in each layers individu- 
ally. 

The radiative transfer module of the original code treated 
scattering on hydrocarbon hazes in the atmosphere (Domagal- 
Goldman et al., 2008; Haqq-Misra et al., 2008). We pro- 
duced new optical data tables for liquid and icy water 
droplets. The new values were calculated for 80 differ- 
ent log-spaced droplet radius values between 0.1 and 100 
microns. 

3. The cloud model 

3.1. Optical properties of cloud droplets 

We assume that the liquid droplets are spherical and 
deformation due to aerodynamical flows around the droplets 
can be neglected. Ices have a crystal structure that de- 
pends on the environment they form in, and the struc- 
ture influences their optical properties (Liou, 1993). These 
shape effects are crucial in Earth atmospheric sciences 
where observational data in good quality is abundant, but 
such effects are secondary in exoplanet atmosphere and 
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paleo-climate problems. Therefore we neglect these shape 
effects and use the equivalent sphere assumption. 

The complex refractive index value of liquid water (Segel- 
stein, 1981) and water ice (Warren and Brandt, 2008) is 
used to calculate the optical properties of droplets as a 
function of radius and wavelength using Mie calculations 
(Bohren and Huffman, 1983). The absorption (Q a bs), scat- 
tering coefficients (Q sca ), and the asymmetry parameter 
(ffsca) of the droplets are shown on Figs. 1, 2, and 3, re- 
spectively. These coefficients are converted into the cross- 
section area of a single droplet: 

a = Q(27rr), (2) 

where a is the (absorption or scattering) cross-section, Q 
is the corresponding coefficient, and r is the droplet radius. 
The asymmetry parameter describes whether the droplet 
is an isotrop scatterer (g sca = 0; light is scattered to all 
directions uniformly), forward scatterer (g sca = 1; analog 
of a perfect lens) or backward scatterer {g sca = — 1; analog 
of a perfect mirror). All of the above describe quantities 
are angle independent for spherical particles. 

The diagonals of Figs. 1, and 2 show that the absorp- 
tion and scattering coefficients are larger than unity, if the 
photon's wavelength is similar to the droplet's radius. In 
this case, the cross-section area of the droplet is several 
times its geometrical area due to a resonance effect be- 
tween the light and the droplet. 

Droplets between 0.2 and 3 microns are strong scatter- 
ers but hardly absorb light. This wavelength range coin- 
cides with the peak of the solar spectral energy distribution 
(SED). However, if the star's SED is more pronounced at 
wavelengths shorter than 0.2 microns, or longer than 3 mi- 
crons, the water droplets become more efficient absorbers. 
The amount of light scattered or absorbed by the clouds 
thus depends on the spectral type of the star as well as on 
the size distribution of droplets. 

The most significant difference between liquid water 
and water ice particles is that ice particles are strong scat- 
terers, but less efficient absorbers in the thermal IR for 
sizes around 10 micron (see Fig. lb). Due to this differ- 
ence, ice particles scatter the outgoing thermal radiation 
more efficiently than liquid droplets. 

3.2. Definitions and basic considerations 

We adopt and use terminology from the meteorology 
community in the next sections, therefore we clarify some 
definitions and explain some concepts. The saturation wa- 
ter vapor pressure is described by the Clausius-Clapeyron 
equation and is a function of temperature. The pressure 
dependence of the saturation vapor pressure can be ne- 
glected under atmospheric conditions. Relative humidity 
is the ratio of partial water vapor pressure to the satura- 
tion vapor pressure. If the relative humidity is larger than 
unity, the air is in a super-saturated state. The amount 
of excess water vapor pressure relative to the saturation 
level is called the level of super-saturation. The saturation 



pressure relative to liquid water and water ice is calculated 
according to Rogers and Yau (1989). If the atmosphere is 
saturated relative liquid water, it is super-saturated rela- 
tive to water ice. 

The process of condensation (or nucleation) can be het- 
erogenous or homogeneous (Rogers and Yau, 1989). Ho- 
mogenous nucleation means that the water vapor molecules 
collide with each other and spontaneously form stable wa- 
ter droplets. This process is not observed in the atmo- 
sphere for two reasons. 1) Homogenous nucleation requires 
extreme levels of super-saturation (several times the satu- 
ration vapor pressure), however the super-saturation levels 
observed in the atmosphere are only a few percentage at 
maximum. 2) Small, sub-micron sized aerosols are always 
present in the atmosphere. The water molecules prefer- 
entially collide with these aerosols and are deposited on 
their surfaces. This process is called heterogenous nucle- 
ation. Therefore, heterogenous nucleation removes the ex- 
cess water vapor before homogenous nucleation could be 
initiated. 

These aerosols or condensation nuclei can be naturally 
produced by erupting volcanoes (ash particles), bursting 
bubbles of salty water (as a bubble on the ocean or sea 
bursts, the tiny droplets are transported away from the 
surface and dried by the wind to form tiny salt crystals), 
dust particles from the land, the agglomeration of various 
sulfate, nitrate, ammonium or organic chemicals, or bac- 
teria and spores can serve as condensation nuclei (Seinfeld 
and Pandis, 2006). Aerosols are also produced by anthro- 
pogenic sources on Earth. The number density and size 
distribution of these aerosols can be measured on Earth 
(see e.g., Meszaros and Vissy, 1974; Hoppel et al., 1990), 
but their properties are unknown on exoplanets. There- 
fore we assume that some amount of abioticly produced 
aerosols are present in all atmospheres such that heteroge- 
nous nucleation can be initiated. 

The formation of cloud droplets can be initiated by 
increasing the relative humidity above the saturation level 
in two ways. 1) Adding moisture to the air, or radiatively 
cooling it without displacing the parcel. This is how fog on 
the ground forms. 2) Increase the altitude of the air parcel. 
In such an environment, the air parcel adiabaticly expands 
and cools until the saturation level is reached after which 
condensation starts. 

Cloud properties on Earth are studied in detail using 
satellite measurements (Rossow et al., 2005; Mace et al., 
2009). Large-scale features of Earth cloud patterns in the 
tropics include the inter-tropical convergence clouds, the 
rising branch of Hadley circulation with typically deep con- 
vective clouds, and the cloud-free descending branch of the 
Hadley circulation. In contrast, mid-latitudes are domi- 
nated by low level clouds. 

Clouds can form as a result of convection during which 
the height of clouds can be larger than their horizontal 
extent (cumulus). Opposed to this, clouds can be hori- 
zontally extended (stratus) or form thin threads (cirrus). 
Stratus and cumulus clouds can be formed at low, and mid- 
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Figure 1: Scattering coefficient of liquid water droplets as a function of droplet radius and wavelength (la), same for water ice (lb). 




Figure 2: Absorption coefficient of liquid water droplets as a function of droplet radius and wavelength (2a), same for water ice (2b). 
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Figure 3: Asymmetry parameter of liquid water droplets as a function of droplet radius and wavelength (3a), same for water ice (3b). 
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altitudes, and they are called alto-stratus or alto-cumulus 
in the latter case. Cirrus clouds are ice clouds thus only 
appear at high altitudes favorable for ice formation. Cir- 
rus clouds can also form via convection (cirrocumulus) or 
form an extended horizontal layer (cirrostratus). There are 
also towering clouds, which occupy more than one height 
ranges (cumulonimbus and nimbostratus clouds). 

There are still many open questions related to cloud 
physics, and we mention only a few of these here. It is an 
active field of research how the relative importance of ra- 
diative, gravity wave, or precipitation processes influence 
convective clouds (Stevens, 2005). It is not fully under- 
stood how the formation and life time of cirrus clouds is 
influenced by different nucleation and radiative processes 
(Jensen et al., 1996; Spichtinger and Cziczo, 2010). Per- 
haps low level stratus clouds are the least well understood 
because their properties (such as the height of the cloud 
base) change by small and slowly varying perturbations in 
radiative forcing, precipitation, and evaporation of falling 
droplets (Dupont et al., 2012). 

It is important to note that although these questions 
are crucial to study and predict the climate on Earth, 
most of these questions represent higher order uncertain- 
ties when dealing with exoplanets and paleo-climate prob- 
lems. The cloud model described in this paper focuses on 
the first order climatic effects of clouds, and it represents 
the first radiative-convective, chemical, and cloud micro- 
physical model developed for exoplanets in the habitable 
zone (HZ). 

3.3. The microphysical cloud model: convective water clouds 

Microphysical cloud models were developed for Earth 
based on observational and laboratory data is available 
that constrain theoretical models (see e.g., Rogers and 
Yau, 1989; Pruppacher and Klett, 1997; Khain et al., 2000; 
Seinfeld and Pandis, 2006). Two issues are important 
when developing a cloud model for exoplanets. 1) Calcu- 
lating the vertical motion, condensation and coalescence 
of a droplet population sclf-consistcntly is computation- 
ally time consuming even in the ID case. 2) Detailed ob- 
servational data is not yet available for super-Earth and 
Earth-like exoplanets to constrain the initial conditions. 

We therefore limit our microphysical model to resolve 
some basic properties of clouds without attempting to sim- 
ulate cloud microphysics fully. These basic properties are 
the following: 

• The cloud base is located at the height where the 
local water mixing ratio equals the saturation mixing 
ratio depending on pressure and temperature. 

• As the convective parcel ascends, it cools and more 
water vapor condenses out. Therefore the size of 
droplets increases with height (Warner, 1969; Rogers 
and Yau, 1989). 



• The cloud top is located at the height where the 
relative velocity between the droplets is large enough 
to initiate coalescence. 

When coalescence becomes efficient, rain drops rapidly 
form and fall towards the surface. We do not directly 
simulate precipitation, only the formation of droplets be- 
cause rain drops are removed from the atmosphere in a 
matter of a few minutes owing to their large terminal ve- 
locities (Rossow, 1978). Therefore the radiative effect of 
rain drops is negligible as they only spend a short amount 
of time in the atmosphere compared to the cloud droplets, 
which remain aloft. 

The falling rain drops however influence the number 
density of cloud droplets as rain drops sweep up the cloud 
droplets as they fall through the cloud (Pinsky et al., 
1999). Thus the number density of cloud droplets is re- 
duced. This effect is considered in our model by the use 
of the precipitation efficiency parameter. 

We consider two vertical columns. In one column, the 
air moves only upwards and clouds form by condensation. 
The other column undergoes downward motion, therefore 
it is expected to be cloud free. 

The initial conditions for the cloudy column are the 
pressure-temperature structure and atmospheric compo- 
sition in the convective zone (T, and Pi, where i is the 
index of layers, i — being on the surface), relative hu- 
midity (x S urf or xo) an d the number density of aerosols 
(^surf) on the surface. From these quantities, we calcu- 
late the height of the cloud deck and top, the radius and 
number density of droplets as a function of height. We 
assume single-sized droplets at a given height, thus we do 
not follow the evolution of a droplet size distribution. This 
simplification is justified because the effective optical prop- 
erties of a droplet size distribution is well represented by 
a single-sized droplet having the equivalent radius of the 
distribution (Hu and Stamnes, 1993). 

An air parcel on the surface becomes warmer than its 
environment, therefore it starts to rise and adiabaticly ex- 
pands. If we assume that the air parcel does not exchange 
material with the surrounding air, the mixing ratio of wa- 
ter vapor will be constant as a function of height (see how- 
ever Houzc, 1993, on the importance of entrainment). The 
number density of aerosols on the other hand will decrease 
as the parcel expands: 

Pi 

n i = ^surfTJ i (3) 
-*surf 

where rn is the number density of aerosols at pressure Pi 
above the surface, P SU rf is the pressure at the surface. The 
pressure scale height of the atmosphere is 10 km, therefore 
the number density of droplets only slightly decrease at 
the typical height of low level liquid water clouds, but the 
modification is significant for high cirrus clouds which form 
above 10 km. 

When the mixing ratio of water in the parcel becomes 
larger than the saturation mixing ratio at the given tem- 
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perature (as described by Rogers and Yau, 1989), the ex- 
cess water vapor condenses onto the surface of cloud con- 
densation nuclei (CCN). The first layer where this happens 
marks the deck of the cloud layer. Above the cloud deck, 
the mixing ratio of water is equal to the saturation mixing 
ratio in that layer. The excess water vapor density in layer 
i {pex,i) is the difference between the water vapor density 
in layer i — 1 and the saturation vapor density in layer 
i. The vapor density is distributed uniformly amongst the 
droplets in the grid cell, and the mass of the water droplets 
in layer i is 

m i =m i _i+ Pex,i/ n i- ( 4 ) 

Initially (below the cloud deck), the droplet mass is zero. 
The radius of the droplets as a function of height can be 
calculated and it increases with height in agreement with 
observations. 

The Reynolds number of the droplets is 



where r, is the droplet radius, p is the density of water (1 
g/cm -3 ), vt,i is the terminal (or fall) velocity of the droplet 
(Rossow, 1978), and fa is the dynamic viscosity of the gas 
(Rogers and Yau, 1989). The Reynolds number describes 
the aerodynamical flow around the droplet. If it is small, 
the flow is laminar, high Reynolds number indicates a tur- 
bulent flow around the droplet. Once the flow becomes 
turbulent, the relative velocity between two droplets is in- 
creased and thus coalescence becomes efficient (Rossow, 
1978). The critical Reynolds number above which coa- 
lescence starts is chosen to be 200. Therefore, the cloud 
top is located at the height where the Reynolds number of 
the droplet becomes larger than 200. There are no droplets 
present above this height, as they would efficiently increase 
in radius due to collisions, and precipitate. The value of 
200 was chosen because the droplet sizes at the cloud top 
is similar to the representative droplet sizes observed on 
Earth. We discuss the effects of different critical Reynolds 
numbers in Sect. 6.1. 

Falling rain drops sweep up cloud droplets during their 
descent. An average cloud droplet has a size of = 10 
microns, an average rain drop is rn — 1 mm. The mass 
ratio of a cloud droplet and a rain drop, assuming spher- 
ical particles, is r^/rf-, = 10~ 6 , thus 10 6 cloud droplets 
must merge to produce a single rain drop. The number 
density of cloud droplets is 100 cm~ 3 , and the number 
density of rain drops is 10 -4 (Rogers and Yau, 1989). The 
ratio of number densities is also 10~ 6 , thus the precipi- 
tation efficiency is of the order of unity. As we do not 
follow the coalescence of a droplet distribution, we intro- 
duce a parameter called precipitation efficiency, e p , which 
can vary between zero and unity. The number density of 
cloud droplets is reduced by a factor of (1 — e p ) due to 
precipitating rain drops. 

The free parameters that determine this cloud layer 
are the surface relative humidity and CCN concentration, 



liquid water cloud fraction, and precipitation efficiency. 

This convective cloud layer is self-consistently calcu- 
lated in our ID model and it is physically motivated. Such 
a cloud model provides a useful first order approximation 
on exoplanetary cloud properties, where the observational 
data in not abundant. However, our simulations show that 
for Earth, where detailed observations of cloud properties 
exist, the resulting convective cloud layer for humidity and 
CCNs cannot reproduce both the average surface temper- 
ature and global albedo of Earth. We find that a precipi- 
tation efficiency of 0.8 is necessary to obtain 289 K surface 
temperature, but the planetary albedo is only 0.21. A pre- 
cipitation efficiency of 0.2 is needed to obtain a planetary 
albedo of 0.3, but the surface temperature is 275 K then. 
The reason why the surface temperature and albedo are 
not reproduced with one cloud layer is because the cooling 
effect of the low level water cloud dominates. If the single 
cloud layer were to form higher in the atmosphere, both 
the albedo and surface temperature could be reproduced. 
However, the only handle we have on the height of the 
convective cloud layer is through the surface relative hu- 
midity. Thus, in order to raise the height of the convective 
layer, we need to reduce the surface relative humidity and 
adopt a value which is different from what is observed on 
Earth. 

In order to reproduce both the globally averaged sur- 
face albedo, temperature and the energy budget of Earth, 
we introduce a second cloud layer: cirrus clouds which 
provide additional warming effect. 

3.4- Cirrus clouds on Earth 

The cloud deck of cirrus clouds cannot be determined 
self-consistently in a ID model because it is a result of 
horizontal motions (see Sect. 3.2). Therefore we use a 
paramctrization to determine the height of the cloud layer. 
Cirrus clouds on Earth form typically at temperatures of 
230 K, which corresponds to a height of 12 km (Liou, 1986; 
Dowling, 1990). There are two possible ways to proceed. 
One can either scale the cirrus height relative to the height 
of the troposphere. The other possibility is to place the 
cloud deck where the temperature is 230 K. Clouds do 
not form above the troposphere, thus the advantage of the 
former method is that it expresses the cloud height as a 
function of its natural maximum value. The advantage of 
the latter method is that a fixed cloud temperature cor- 
responds to a fixed super saturation relative to water ice, 
and it also directly informs us about the greenhouse effect 
of these clouds. Neither of the two methods have an over- 
whelming advantage over the other. However, we found 
in our simulations that linking the cloud height with the 
tropopause height sometimes results in non-convergence 
and continuous surface cooling. Therefore, for pure nu- 
merical reasons, we use the latter method to determine 
the height of the cloud deck. 

The size of the particles as a function of height is de- 
termined in the next step. We assume that the air parcel 
is saturated relative to liquid water, thus super-saturated 



6 



relative to ice. This assumption is commonly used in cir- 
rus cloud models (Starr, 1985; Starr and Cox, 1985; Liou, 
2000). The difference condenses onto ice condensation nu- 
clei (ICN) and the particle sizes as a function of height can 
be determined. The number density of ICN is calculated 
from the same surface aerosol number density as for liquid 
water clouds for simplicity. 

The height of the cloud top is reached where the parti- 
cle Reynolds number is larger than 200, as for liquid water 
droplets. The falling particles sweep up other particles 
as they settle thus reducing the particle number density. 
Again, for simplicity we assume the same precipitation ef- 
ficiency as for liquid water clouds. 

As a result of our assumptions, we only need one ad- 
ditional free parameter (the ice cloud fraction) to describe 
the second cloud layer. 

Once the particle properties in the two cloud layers are 
determined, we calculate the visible and IR fluxes. We 
use flux averaging to determine the equilibrium pressure- 
temperature structure for a cloudy atmosphere (Kitzmann 
et al., 2010; Goldblatt and Zahnle, 2011). The initial or 
previously determined T-P profile is used to calculate the 
radiative fluxes. The T-P profile and mixing ratios are 
kept fixed and the incoming and outgoing fluxes are cal- 
culated for four cases: cloud- free T-P profile (F c f), T-P 
profile with liquid water clouds only (Fi c ), with ice clouds 
only (F ic ), and with both clouds overlapping (F a ). These 
four fluxes are averaged using the cloud fractions as weight 
and assuming a random overlap between the two cloud lay- 
ers: 

F = F cf x (1 - /, - fi + fifi) + F lc x (/, - fifi) + 
F ic x (fi-fifi)+F x (fift), 

where fi is the cloud fraction of the liquid water cloud, 
and fi is the cloud fraction of the ice cloud. The average 
flux (F) is used in Eq. 1 to calculate the new, updated 
T-P profile. 

4. Model set up and tests 

4-1- Initial conditions 

Unless otherwise stated, we use the following initial 
conditions for all simulations performed in this paper. We 
adopt the physical and orbital parameters of the Earth- 
Sun system and initiate the simulations with the stan- 
dard Earth atmosphere composition: = 10~ 2 , fcn 4 = 
1.6 x 10~ 6 , fco 2 = 3.55 x 10~ 4 , fo 2 = 0.21, where / is the 
volume mixing ratio of argon, methane, carbon-dioxide, 
and molecular oxygen respectively. The rest of the atmo- 
sphere consists of molecular nitrogen (/jv 2 = 0.77). Other 
important tracer species are water vapor, and ozone whose 
initial mixing ratio is set to fH 2 o = 10~ 5 in the strato- 
sphere, and fo 3 = 10~ 6 . The relative humidity is a free 
parameter in the code. However, we set it to 0.77 on the 
surface (global average Earth value) for the convergence 
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Figure 4: Convergence of the surface temperature, stratosphere tem- 
perature and cloud heights. The x axis shows the number of iter- 
ations performed by the climate code. The y axis (left side) is the 
temperature. The upper solid line shows the surface temperature, 
the lower solid line shows the stratosphere temperature at 10~ 3 bars, 
where the stratosphere temperature is the highest. The y axis (right 
side) is height in km, where the lower pair of dotted lines show the 
height of the liquid water cloud. The upper pair of dotted lines 
show the ice cloud heights. Convergence in all of these quantities is 
reached after 5 X 10 4 climate iterations. 



tests described in Sect. 4.2. Most other chemical species 
are initiated with / ot hor — 10~ 20 . 

The initial temperature and pressure profile is calcu- 
lated as follows. The surface pressure is set to 1 bar, the 
pressure on the top of the atmosphere is 5 x 10 -5 bar rep- 
resenting about 70 km in current Earth's atmosphere. The 
pressure grid is initialized such that it is uniformly spaced 
(6) m log-pressure coordinates (but see Par. 3 in Sect. 4.2). 
The temperature on the surface is set to 289 K, the tem- 
perature in the stratosphere is initially 200 K. The height 
of the troposphere is set to 10 km (but evolves during the 
simulation) and the temperature in the troposphere drops 
linearly to reach 200 K at the top of the troposphere. Once 
the pressure and temperature is known, the height of the 
layers can be calculated assuming hydrostatic equilibrium 
and using the ideal gas equation of state. 

As discussed in Sect. 2.2, the surface albedo of the 
planet was initially set to a su rf = 0.22 in the cloud-free 
model to mimic the effects of clouds. However, the clear 
sky albedo of Earth is 0.15 (Chapter 5 of Pierrehumbert, 
2011, based on ERBE measurements). The clear sky albedo 
is the combined albedo of the gaseous atmosphere (reflect- 
ing light via Rayleigh scattering) and the reflectance of the 
surface, no clouds arc considered. The clear sky albedo in 
the model can be measured as the ratio of the outgoing 
and incoming visible/near-IR flux on the top of the at- 
mosphere (F^-top / F^top) ■ We set the surface albedo such 
that F-ptop/Fj.top = "dear = 0.15. If we follow this pro- 
cedure, we find that the surface albedo has to be set to 
"surf = 0.13 (in agreement with Trenberth et al., 2009). 
This means that 2% of the incoming radiation is reflected 
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off the atmosphere by atoms and molecufes in the clear sky 
case. The surface temperature is then between 292 and 295 
K depending on the number of vertical layers considered 
(see Sect. 4.2). As this value is larger than the global 
average of 289 K, it already suggests that cooling effect of 
water clouds dominates over the warming effect of cirrus 
clouds on Earth. The zenith angle is set to 60° in these 
simulations. 

The time step in the climate model is adaptive and it 
is regulated by the maximum temperature difference in all 
layers (AT). It is increased by a factor of 1.1, if AT is 
smaller than 5 x 10 ~ 2 K, the time step is divided by two, 
if AT is greater than 10 -1 K. These values are an order of 
magnitude smaller than in the clear sky model due to the 
finer vertical resolution in our model with clouds. 

Clear sky calculations call the climate code 200 times, 
the photochemistry code 400 times, and the whole pro- 
cedure is repeated 2-3 times until convergence is reached 
(Scgura et al., 2010). The cloudy sky model calls the cli- 
mate code 4000 times, the photochemistry 800 times, and 
repeats the procedure 20 times. The significant increase 
in iteration number is necessary because of the small time 
steps and different initial conditions. We illustrate the 
convergence of surface temperature, stratosphere temper- 
ature, and cloud heights on Fig. 4. The cloud model 
parameters of the simulation are described in the third 
convergence test of Sect. 4.2. We see that both the sur- 
face temperature and cloud heights reach convergence af- 
ter 5 x 10 4 climate iterations. However, in some cases 
where the T-P profile differs significantly from the initial 
T-P profile, convergence is reached later, thus iterating 10 5 
times is justified. 

4-2. Convergence tests 

Clear-sky test. Using the above mentioned initial condi- 
tions, we performed clear-sky simulations (no clouds) with 
increasing number of grid cells to verify the convergence 
of the code. The number of layers varies between n z = 60 
to 600 with a step size of 30. Figure 5a shows the surface 
temperature and clear sky albedo as a function of layer 
numbers. As expected, both the surface temperature and 
the clear-sky albedo converges as the number of vertical 
layers increases. The surface temperature difference be- 
tween n z = 570 and 600 runs is 0.01 K, the surface albe- 
dos are basically equal (a difference occurs only at the fifth 
decimal point). Therefore we conclude that convergence is 
reached for n z — 600. 

Figure 5b shows the temperature, ozone and water mix- 
ing ratios as a function of height after convergence was 
reached for n z = 600. The dotted lines show the US Stan- 
dard Atmosphere (1976) profiles. The agreement between 
the standard atmosphere profiles and the derived profiles 
are not perfect in Fig. 5 for two reasons. 1) The standard 
atmosphere profiles are based on globally averaged mea- 
surements, therefore these profiles might not represent an 
atmosphere that is in radiative and chemical equilibrium. 
2) The correct Earth surface temperature is influenced by 



clouds, which is not accounted for in this model atmo- 
sphere. 

Parameterized cloud model. We adopt a homogenous cloud 
model similar to the one used by Goldblatt and Zahnlc 
(2011) for their low level water clouds. We assume a 40% 
cloud coverage, cloud deck at 0.85 bar, cloud top at 0.7 bar, 
mono-disperse droplets meaning that all droplets have the 
same radius, 11 micron, and column density of cloud water 
(or water path) 40 g/m 2 . The number of layers are varied 
between 150 and 600 with a step of 30. We found that 
using less than 150 layers introduce larger uncertainties in 
the surface temperature and albedo because there are only 
one or two cloudy layers, which can potentially extend to 
a much wider or narrower pressure range than 0.85 and 
0.7 bar. The number of cloudy layers are three and twelve 
for 150 and 600 vertical layers, respectively. 

The surface temperature and albedo are shown on Fig. 
6a, the temperature, ozone and water mixing ratios are 
shown on Fig. 6b. The surface temperature is between 
273-285 K, the albedo of the planet is between 0.25-0.27 
depending on the number of layers used. Convergence is 
observed for the surface temperature and for the planetary 
albedo. The cloud deck and top are not exactly at the 
prescribed pressures due to the vertical discretization of 
the pressure grid. As the exact pressure values of the cloud 
deck and top fluctuates, so does the surface temperature 
and albedo values as seen on Fig. 6. 

The microphysical cloud model. In this convergence test, 
we use the microphysical cloud model as described in Sects. 
3.3 and 3.4. The free parameters of the microphysical 
cloud model are the number density of aerosols (n sur f), the 
cloud fractions (/;, and fi), the relative humidity on the 
surface (xsm-f), and the precipitation efficiency (e p ). Here 
we use n sur f = 100 cm~ 3 (Miles et al., 2000), Xsurf = 0.77, 
which both are average measured Earth quantities. Fi- 
nally, e p — 0.8, /; = 0.4, and fi — 0.25 are used to obtain 
surface temperatures equal to the measured average Earth 
value of 289 K. 

In a pressure grid which is uniformly spaced in log pres- 
sure (see Sect. 4.1), there are only two cloudy layers even 
for 600 grid cells, thus we found no real convergence. One 
could increase the number of layers above 600, but that 
renders the computations too expensive. We change the 
pressure grid instead. The ratio of grid spacing at the top 
to the grid spacing at the bottom of the atmosphere is 
set to 15, which means that the lower atmosphere (where 
clouds are located) is well resolved, the grid spacing at the 
top of the atmosphere is coarse. Under such conditions, 
there are on average 4 and 10 cloudy grid cell per cloud 
layer for 150 and 600 grid cells, respectively (see Fig. 7b). 

Figure 7a shows the surface temperature and albedo 
as a function of layer number. The surface temperature 
varies between 277 and 289 K, while the albedo varies be- 
tween 0.345 and 0.3. The albedo decreases and surface 
temperature increases with greater layer numbers. Figure 
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Figure 5: Left side: convergence test for the clear-sky case. The solid line shows the surface temperature as a function of the number of layers 
used in the model. The dotted line shows the surface albedo. Right side: the temperature, ozone and water mixing ratios as a function of 
height for n z = 600. The dotted lines show the US Standard Atmosphere (1976) profiles. 




Figure 6: Left side: the surface temperature and albedo as a function of layer number used for the parameterized cloud convergence test. 
We assume a 40% cloud coverage, the cloud deck and top are located at P = 0.85 and 0.7 bars, respectively, the droplets have a radius of 
11 microns, and the total water path is 40 g/m 2 (Goldblatt and Zahnlc, 2011). Right side: the temperature, water and ozone mixing ratio 
profiles as a function of height for n z = 600. The dotted lines show the US Standard Atmosphere (1976) profiles. 
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Figure 7: Left side: the surface temperature and albedo as a function of layer number used for the convergence test of the microphysical 
cloud model. We assume a surface relative humidity of 77%, an aerosol number density of 100 cm~ 3 , a liquid water cloud fraction of 40%, 
cirrus cloud fraction of 25%, and a precipitation efficiency of 0.8. Right side: the properties of the cloud layers. The y axis shows the height 
of the layer, the color represents the radius of the droplets in the layers. The red contours show the liquid water cloud properties, the blue 
contours correspond to the ice cloud. 



7b shows the properties of the clouds in the simulations. 
The y axis shows the height of the cloudy layers, the cen- 
ter of the layers is marked with a '+' sign, and the col- 
ors correspond to the droplet radius. Red contours show 
the properties of the liquid water cloud, the blue contours 
show the ice cloud properties. In agreement with observa- 
tion and theory, the droplet size is increasing with height 
(e.g., Warner, 1969; Rogers and Yau, 1989). There are no 
droplets with a radius above 20 microns at the adopted 
pressure and temperature conditions because of rain-out. 

We use 600 layers and a grid spacing ratio of 15 for the 
following simulations. 

4-3. CO2 doubling 

We perform one simulation to measure the climate sen- 
sitivity due to CO2 doubling. This is another way to vali- 
date our model. We use the standard cloud parameters for 
Earth (see Sect. 5.1 with precipitation efficiency 0.8) but 
with a doubled CO2 concentration (fco 2 = 7.1 x 10~ 4 ). 
We find that the surface temperature is 292.5 K, and the 
planetary albedo is 0.29. The present Earth values are 289 
K and 0.3. Both cloud layers form at slightly increased al- 
titudes. The cloud decks are located at 1.6 km and 13.1 
km. These values in our present Earth simulation are 1.3 
km and 12.1 km. 

The climate sensitivity of our model to CO2 doubling is 
thus 3.5 K which compares favorably to the IPCC best es- 
timate of 3 K with a likely range of 2 - 4.5 K (IPCC, 2007). 
Therefore we conclude that our model reliably reproduces 
the climate of Earth. 

5. Results 

We perform five sets of simulations and in each of these 
sets, we keep the value of four parameters constant and 



vary only the fifth one. We could map the whole param- 
eter space by simulating all possible combinations of the 
free parameters, but we feel that our understanding is bet- 
ter served by varying only one parameter at a time and 
exploring the effects individually. 

5.1. Precipitation efficiency 

First, we set the number density of aerosols to 100 
cm" 3 , the relative humidity to 77%, the liquid water cloud 
fraction to 0.4, and the ice cloud fraction to 0.25 in agree- 
ment with Goldblatt and Zahnle (2011), and vary the pre- 
cipitation efficiency between and 1 with a step of 0.05. 
This set of simulations is used to determine what precipi- 
tation efficiency is needed to reproduce the average surface 
temperature, albedo and GEB of Earth by adopting the 
measured values of the other four cloud parameters. 

We show the albedo and the surface temperature as a 
function of precipitation efficiency on Fig. 8a. A low pre- 
cipitation efficiency produces geometrically and optically 
thick clouds and the albedo effect of clouds dominates over 
their greenhouse effect. As the precipitation efficiency is 
increased, the clouds become more optically thin, more of 
the incoming solar radiation reaches and warms the sur- 
face. For a precipitation efficiency of unity, there are no 
clouds present and we obtain the clear sky albedo and sur- 
face temperature values. The precipitation efficiency value 
resulting in 289 K surface temperature is 0.8, and the cor- 
responding albedo is 0.3 in agreement with observations 
(Chapter 5 of Pierrehumbert, 2011, based on ERBE mea- 
surements) . 

The cloud properties are illustrated on Fig. 8b, which 
shows the height of the cloudy grids as a function of precip- 
itation efficiency with '+' signs, and the colors represent 
the droplet radius in microns separately for the liquid and 
ice clouds. As the surface temperature increases, so does 
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Figure 8: The precipitation efficiency is varied between and 1 for /; = 40%, fi = 25%, aerosol number density of 100 cm -3 , and relative 
humidity of 77%. Left side: the albedo and surface temperature as a function of precipitation efficiency. Right side: the height of the cloudy 
layers with the '+' signs. The color indicates the droplet sizes. 



Table 1: Global energy budget. The first column describes the flux 
type, the second column is the observed flux values from Zhang et al. 
(2004), the third column is the flux values calculated in our model 
to reproduce Earth's GEB. The relative error between the measured 
and simulated quantities is a few percent at maximum, therefore we 
conclude that the agreement is good. 



Measured Simulation 



[W/m 2 ] [W/m 2 ] 

Incoming solar flux 341.8 340.3 

Absorbed by atoms 70.9 73.1 

Reflected by atoms and clouds 81.7 76.6 

Reached the surface 189.2 190.6 

Leaves the surface 24.0 24.8 

Leaves the atmosphere 105.7 101.4 

Emitted IR flux from the surface 395.6 396.8 

Back radiation to the surface 344.7 344.8 

IR flux leaving the atmosphere 233.3 238.8 



the height of the cloudy grid cells. The size of the droplets 
increases with height in accordance with observations and 
theory and their maximum size does not exceed 20 microns 
due to the onset of precipitation. 

The GEB describes two processes: 1) how the incom- 
ing solar flux is absorbed by atoms, reflected by clouds, 
how much solar flux reaches the surface, leaves the surface 
and finally the atmosphere; 2) the amount of thermal IR 
flux emitted by the surface, how much IR flux is back radi- 
ated to the surface, and leaves the atmosphere. The GEB 
of Earth is determined by satellite observations (Zhang 
et al., 2004) and we compare the GEB of our simulation 
to the observed values in Tab. 1 for a precipitation effi- 
ciency of 0.8. The maximum relative error of the observed 
and measured quantities is 6%, the typical relative error is 
3%, thus our model reproduces Earth's GEB with a pre- 
cipitation efficiency of 0.8. 



5.2. Cloud fractions 

Clouds have two competing climatic effects. They re- 
flect the incoming solar radiation thus cool the surface, 
and they scatter as well as absorb and re-emit the outgo- 
ing IR radiation. Whether the net climatic effect of a cloud 
layer is positive or negative depends on the droplet mate- 
rial and size, as well as the temperature at the cloud's 
location. If the cloud layer is located at high altitudes, 
thus the cloud temperature is low compared to the surface 
temperature, the cloud has a net warming effect on the 
surface. Low altitude clouds with temperatures similar to 
the surface temperature tend to cool the surface. These 
considerations are verified by our model. 

We use the measured aerosols number density n sur f = 
100 cm~ 3 , relative humidity (xsurf = 0.77) of Earth, set 
the precipitation efficiency to e p — 0.8 as derived in the 
previous section, set fi to be 0.25, and vary the liquid 
water cloud fraction (/;) between and 1 with a step of 
0.05. 

Figure 9a shows the albedo and surface temperature 
as a function of cloud fraction. Even when the liquid 
water cloud fraction is zero, there is an ice cloud layer 
present. As expected, the albedo increases (surface tem- 
perature decreases) with larger cloud fraction. Figure 9a 
shows a linear correlation between the albedo and cloud 
fraction up-to 90% cloud cover. The albedo and surface 
temperature for 100% cloud coverage are 0.57 and 243 K, 
respectively. We also see that the surface temperature and 
albedo levels off above a cloud fraction of 0.9. This behav- 
ior is explained in Fig. 9b, which shows the cloud prop- 
erties, height of the cloudy layers and droplet radius. As 
the surface temperature decreases, the height of especially 
the ice cloud is reduced. As the temperature difference 
between the surface and at the location of the ice cloud 
(230 K) decreases, the warming effect of the ice layer van- 
ishes, thus the surface temperature is drastically reduced. 
Above a cloud fraction of 90%, the two cloud layers form 
one single layer and the albedo is only slightly raised by 
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Figure 9: The liquid water cloud fraction is varied between and 100% for /j = 25%, an aerosol number density of 100 cm -3 , relative 
humidity of 77%, and precipitation efficiency of 0.8. Left side: the albedo and the surface temperature as a function of cloud fraction. Right 
side: the height of the cloudy layers and the size of the droplets. 
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Figure 10: The water ice cloud fraction is varied between and 100% for /; = 40%, an aerosol number density of 100 cm -3 , relative humidity 
of 77%, and precipitation efficiency of 0.8. Left side: the albedo and the surface temperature as a function of cloud fraction. Right side: the 
height of the cloudy layers and the size of the droplets. 
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Figure 11: The number density of aerosols is varied between 70 cm 3 and 300 cm 3 for /; = 40%, /; = 25%, relative humidity of 77%, and 
precipitation efficiency of 0.8. The albedo and surface temperature are shown on the left, and the cloud properties on the right. 
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Figure 12: The relative humidity is varied between and 100% for a liquid water cloud fraction of 40%, cirrus cloud fraction of 25%, aerosol 
number density of 100 cm -3 , and precipitation efficiency of 0.8. The albedo and surface temperature are shown on the left, and the cloud 
properties on the right. The y axis range differs from Figs. 9b and lib. The dotted lines indicate the range values used on those figures. 



further increasing the cloud fraction. 

In the following, we set /; to 0.4 and vary the cirrus 
cloud fraction between and 1 with a step of 0.05. Ini- 
tially, the surface temperature increases with higher cloud 
fraction and warms the surface (see Fig. 10a). As a re- 
sult, cirrus clouds form at larger altitudes, as their height 
is fixed at 230 K in this model (see Fig. 10b). When clouds 
form at lower pressures, the number density of ICNs and 
cloud particles decrease in accordance to Eq. 3 and the 
clouds become optically thinner. As a result, the warming 
effect of cirrus clouds levels off above fi > 0.5. 

As seen on Fig. 10b, cirrus clouds form at a height of 
10-20 km. The pressure at 20 km is 7 x 10~ 2 bar, thus 
the initial number density of aerosols at this height is 7 
cm" 3 according to Eq. 3. This value is further reduced to 
1.4 cm -3 for a precipitation efficiency of 0.8. The number 
density of particles for a cirrus cloud at a height of 10 km 
is roughly four times higher. Thus cirrus clouds at a height 
of 20 km are thin and have a low greenhouse effect. 

5.3. Number density of aerosols 

In this section, we set the fi to 0.4, fi to 0.25, the 
relative humidity to 0.77, the precipitation efficiency to 
0.8, and change the number density of aerosols between 70 
cm" 3 and 300 cm" 3 with a step of 10 cm" 3 . The number 
density of 70 cm -3 is measured at oceanic areas of Earth, 
while a number density of 300 cm" 3 is measured in heavily 
populated industrial areas (Miles et al., 2000). 

Figure 11a shows that the albedo and the surface tem- 
perature are also sensitive to the aerosol number density. 
If the number density of aerosols was zero, there would 
be no clouds on the planet, and the surface temperature 
would be 295 K with an albedo of 0.15, the clear sky case. 
We find that high aerosol number density produces geo- 
metrically (and optically) thick clouds, thus increasing the 
albedo and reducing the surface temperature of the planet. 
If the number density of aerosols is high, the same amount 



of water vapor can condense onto more aerosols, thus the 
radius of cloud droplets gradually rises with height and the 
cloud layer is thick (see Fig. lib). If the number density 
of aerosols is low, the cloud layer becomes thin and the 
droplet size steeply increase with height. 

5.4- Relative humidity 

We set the number density of aerosols to 100 cm -3 , the 
fi to 0.4, fi to 0.25, the precipitation efficiency to 0.8, and 
vary the surface relative humidity between and 1 with a 
step of 0.05. This way we consider the full range between 
desert and ocean planets. 

Figure 12a shows that a relative humidity of 0% cools 
the planet to 273 K. In this case, no clouds are present in 
the atmosphere. The surface temperature drops because 
an efficient greenhouse gas, water vapor, is removed from 
the atmosphere. The simulation with 5%, 10% and 15% 
relative humidity are excluded because the cumulus clouds 
form at altitudes larger than the cirrus clouds, which is 
unphysical. Therefore, our cloud model is not applicable 
for low relative humidities. 

The climate is insensitive to changes in the relative hu- 
midity compared to the other cloud parameters. As the 
convective water cloud is located where an ascending air 
parcel becomes saturated with water vapor, higher rela- 
tive humidity decreases this height. If the relative hu- 
midity is small, clouds form at high altitudes (see Fig. 
12b) and their green house effect dominates over their 
albedo effect compensating for the lack of water vapor 
in the atmosphere. At high relative humidities, the con- 
vective clouds form close to the surface and thus have a 
reduced greenhouse effect. However, the atmosphere con- 
tains large amount of water vapor to compensate. If the 
relative humidity is 100%, convective clouds form already 
at ground level, as any vertical displacement results in 
super-saturation. 
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6. Discussion 

6.1. Critical view on the cloud model 

Often convection is treated as a ID diffusion process in 
cloud models for brown dwarfs and giant planets (see e.g., 
Ackerman and Mar ley, 2001). The formed cloud droplets 
can be displaced to both directions and this process ar- 
tificially increases the vertical extent of the cloud layer. 
On Earth, we see that the regions where convection is 
upward and downward are physically separated, thus the 
downward motion does not influence the cloud properties. 
Therefore we adopted this concept in our microphysical 
model. 

Cirrus clouds are produced by overrunning and under- 
cutting on Earth (Sect. 3.2). However, such horizontal 
motions cannot be self-consistcntly modeled in ID codes, 
thus we are left to parameterize the height of the cloud 
deck. The adopted parametrization is probably the largest 
uncertainty in our model because it is solely based on 
Earth observations. The temperature or super-saturation 
level where cirrus clouds form could change, if the rela- 
tive humidity, the circulation pattern of the atmosphere, 
etc changes. Unfortunately we are not aware of any cirrus 
model which could take such effects into account in a ID 
model. 

We consider cloud formation only by condensation, and 
consider one representative droplet size at a given height 
because the optical properties of a droplet size distribu- 
tion can be approximated by a single droplet size (Hu and 
Stamnes, 1993). We do not simulate coalescence and the 
vertical motion of a droplet population. The effects of 
such microphysical processes are included in a simplified 
way in the precipitation efficiency. This parameter could 
in principle be constrained self-consistently by simulating 
the vertical motion, condensation, and coalescence of a 
droplet population. However, such calculation would in- 
troduce more free parameters (e.g., the distribution and 
material of aerosols) and add a significant computational 
overhead to our numerical methods. To keep the free pa- 
rameters limited and the computational time short, we 
introduce the precipitation efficiency. 

The number density of aerosols is determined by vari- 
ous processes such as bubble bursts on the sea, wind lifting 
up dust from the continents, volcanos, etc. (see Sect. 3.2). 
The total number of cloudy layers and the cloud fraction 
per cloud layer are determined by the circulation of the 
atmosphere. For example the rotation period, the friction 
between the surface and air, topography, relative humid- 
ity, number density of aerosols all influence the cloud frac- 
tion. The relative humidity is mostly determined by the 
topography of the planet, by the zonally averaged land/sea 
distribution and surface temperature profile. Detailed two 
dimensional cloud resolving models show that the precip- 
itation efficiency depends on the water vapor content of 
the atmosphere, surface evaporation rate, amount of pre- 
cipitating particles, zonal and vertical wind speeds, etc. 
(Sui et al., 2007; Gao and Li, 2011). As we either do not 



know the above described planet properties or cannot sim- 
ulate the physical processes for exoplanets, we varied these 
parameters independently in Sect. 5. GCMs are able to 
resolve some of the correlations mentioned above, but such 
models are currently too computationally expensive to be 
used for parameter exploration. Therefore we were left to 
treat these parameters independently for now. 

We use the critical Reynolds number in both cloud 
layers to limit the width of the cloud layers. A value 
of 200 is chosen for both layers, but other values could 
also be adopted. The critical Reynolds number deter- 
mines the width of the cloud layer. If Re is smaller (larger) 
than 200, the cloud layers become thiner (thicker) and a 
smallcr(larger) value of precipitation efficiency would also 
reproduce the global average temperature and albedo of 
Earth. However, the critical Reynolds number should not 
be smaller than 70, because the flow around the droplets 
is laminar and coalescence is not initiated in that case 
(Rossow, 1978). Adopting too large values of Re makes the 
cloud too thick geometrically and optically, and a precipi- 
tation efficiency of almost unity is necessary to reproduce 
the observed properties of Earth. We verified by simula- 
tions that the results are not sensitive to the exact value 
of Re. Varying the other cloud parameters yield to large 
trends, and the effect of different Re is negligible compared 
to these trends. 

We do not model the ice albedo feedback. The sur- 
face albedo decreases for higher surface temperature as 
the size of ice-covered areas become smaller due to melt- 
ing. If the surface temperature decreases, the size of ice 
covered areas would increase, which in turn would fur- 
ther decrease the surface temperature. This is a positive 
feedback, as some cooling(warming) promotes further cool- 
ing(warming) of the surface. However, we keep the surface 
albedo fixed in all runs. 

All our model runs were done with a fixed composi- 
tion of non-condensable gases and a surface pressure of 
1 bar. Comparison with Venus, Mars, and Titan indi- 
cates that the composition and mass of terrestrial planet 
atmospheres varies widely. Similar is expected for exo- 
planets. For Earth, the composition (Sagan and Mullen, 
1972; Owen et al., 1979) and possibly the mass (Goldblatt 
et al., 2009) of the atmosphere have evolved with the geo- 
chemical evolution of the planet. There are geochemical 
feedbacks between temperature and composition. Most fa- 
mous is the temperature dependence of CO2 concentration 
by silicate weathering, which is thought to exert a negative 
feedback on the global temperature on 10 5 yrs timescales 
(Walker et al., 1981). Until we have high resolution spec- 
tra for terrestrial exoplanets, the atmospheric composition 
will be difficult to constrain - and modeling such effects is 
thus beyond the scope of our climate model. 

6.2. Comparison of microphysical model to parameterized 
cloud models 

Our standard simulation is described in Sect. 5.1 with 
a precipitation efficiency of 0.8. We compare this simula- 
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tion to the standard simulations of Goldblatt and Zahnlc 
(2011) and Kitzmann ct al. (2010). Both of these cloud 
models are parameterized and have been applied to Earth 
as a test case. 

Kitzmann et al. (2010) uses two cloud layers, like our 
model. Their model reproduces the global average tem- 
perature of Earth, but the planetary albedo is somewhat 
smaller (0.27) than the observed value of 0.3. The plane- 
tary albedo is calculated as the ratio of incoming to outgo- 
ing solar flux in the GEB. This implies that the outgoing 
solar flux of Kitzmann et al. (2010) is 10% smaller than 
the observed value. Kitzmann et al. (2010) adopts 52 ver- 
tical grid cells and their cloud layer does not entirely fill 
one grid cell. The thickness of the grid cell containing the 
low level water cloud is for example 0.9 km, however the 
thickness of the cloud layer is 0.15 km (personal commu- 
nication). Kitzmann et al. (2010) adopts a droplet size 
distribution, number density, and cloud optical depth in 
accordance to observations. The cloud thickness is then 
derived from these quantities. The grid cell containing the 
cloud is determined by the observed pressure at the cloud 
top. 

The cloud heights (or pressures) of Goldblatt and Zahnle 
(2011) are derived from observations. They adopt 30 ver- 
tical layers mostly centered around the troposphere and 
use three clouds layers. The low level liquid water cloud 
is resolved by 3 grid cells, the middle layer by 2, and the 
cirrus cloud layer by 1 grid cell. Their GEB values are 
shown in Fig. 5 of their paper and are in good agreement 
with the observed values. The largest relative error is 5% 
in their GEB fluxes with a typical error of 2-3%. 

Kitzmann et al. (2010) and our work adopt an initial 
droplet number density in agreement with observations. 
Our droplet number density is then reduced due to precip- 
itation. The number density of droplets in Goldblatt and 
Zahnlc (2011) is calculated using the liquid water path, 
average droplet size, the height of the cloud deck and top. 
Their lower liquid water cloud has for example 7 droplets 
per cm 3 , which is lower than in the other two models. This 
difference is explained by the large vertical extent of their 
clouds (1.4 km, compared to 0.4 km in our standard model 
and 0.15 km in Kitzmann et al. (2010)). 

Finally, the results of parameterized cloud models dif- 
fer in details from the results of the microphysical cloud 
model but all three models reproduce the GEB of Earth 
much better than clear sky models with an increased sur- 
face albedo. The advantage of our microphysical cloud 
model over parameterized models is that the cloud proper- 
ties are determined self-consistently (except for the cirrus 
height), thus it is applicable to a much wider range of at- 
mospheric and planetary properties than a parameterized 
model, which is tuned to Earth. The cloud model predict 
how cloud properties change and influence the climate, if 
relative humidity, cloud fraction, aerosol number density, 
atmospheric composition, etc is different than what is ob- 
served on Earth. Parameterized models have difficulty to 
cope with such issues. 



6.3. Sensitivity of the climate 

The simulations performed in Sect. 5 map how the 
structure of the atmosphere is influenced by cloud fraction, 
aerosol number density, relative humidity and precipita- 
tion efficiency. For precipitation efficiency between and 
100%, the albedo and surface temperature varies between 
0.47-0.15, and 259 K - 230 K, respectively. For a liquid 
water cloud fraction between and 100%, the albedo and 
surface temperature varies between 0.16 - 0.58, and 305 K 
- 242 K, respectively. For an ice cloud fraction of and 
100%, the albedo and surface temperature varies between 
0.26 - 0.3, and 282 K - 300 K, respectively. For aerosol 
number density between 70 cm 3 and 300 cm~ 3 , the albedo 
and surface temperature changes between 0.25 - 0.47, and 
293 K - 262 K, respectively. The climate is least sensitive 
to the relative humidity. For relative humidity between 
20% and 100%, the albedo and surface temperatures are 
0.3 - 0.27, and 287 K - 294 K, respectively. A cloud-free 
atmosphere with zero relative humidity has a surface tem- 
perature of 273 K, and albedo of 0.17. 

However, the relative humidity, cloud fraction and pre- 
cipitation efficiency are not independent parameters. The 
relative humidity determines the amount of water vapor 
available for condensation and cloud formation. Therefore, 
assuming a zero cloud fraction and 77% relative humidity 
in Sect. 5.2 is probably unrealistic, if CCNs exist on the 
planet. Assuming less than 20% relative humidity and 
40% cloud fraction in Sect. 5.4, and a precipitation effi- 
ciency of both zero and unity in Sect. 5.1 are unlikely. The 
relative humidity and cloud fraction arc dependent param- 
eters, however the dependency is difficult to constrain even 
with global circulation models. Therefore we treat them 
as independent parameters. Determining the precipitation 
efficiency self-consistently would add new free parameters 
and significant computational overhead to the calculation 
(see Sect. 6.1). 

It shows the intricacy of the field that the albedo and 
climate depends rather strongly on the aerosols. It is 
also observed on Earth that the cloud properties are dif- 
ferent above land, where the aerosol number density is 
high, and above sea, where the aerosol number density 
is low (Rossow and Schiffer, 1999). Models to calculate 
the aerosol number density exists for Earth (Seinfeld and 
Pandis, 2006) and it might be beneficial in the future to 
generalize those models for exoplanets. 

7. Summary 

We developed a microphysical cloud model that can 
self-consistently calculate the cloud heights and droplet 
sizes in ID atmosphere models. Our treatment of low level 
water clouds is inspired by vertical convection where warm 
and humid air rises from the surface and forms clouds, 
and dry air descends at some other location. The height 
of water ice clouds (cirrus clouds) is determined by a pa- 
rameterization in temperature. The free parameters of the 
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model are the surface relative humidity, number density of 
cloud condensation nuclei, cloud fractions, precipitation 
efficiency (and the temperature where cirrus clouds form, 
if Earth is considered). 

We apply our model to Earth. As a first step, we 
adopted the observed relative humidity, aerosol number 
density, and cloud fraction values to determine what pre- 
cipitation efficiency is necessary in the framework of our 
model to reproduce the observed average surface temper- 
ature, albedo, and GEB of Earth. We find that a 80% 
precipitation efficiency is necessary to reproduce the mea- 
sured GEB with good accuracy. Thus both the parameter- 
ized cloud models and our cloud model can reproduce the 
globally averaged observed properties of Earth well. The 
Goldblatt and Zahnle (2011) model reproduces the GEB 
with a relative error of 2-3%, our model with 3-4%, and 
the Kitzmann et al. (2010) model has a relative error of 
10% on the visible flux leaving the atmosphere. 

The advantage of our model over the parameterized 
models is that we calculate the cloud properties self-consistentljGoldblatt, C, Zahnle, K. J., Mar. 2011. Clouds and the Faint Young 
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(expect for the height of the cirrus layer), thus we can 
explore how the cloud properties and the climate is in- 
fluenced by varying the cloud fractions, relative humidity, 
aerosol number density, atmospheric composition, etc. A 
parameterized model is tuned for an exact Earth-analog 
and it is challenging to extend it to a wider parameter 
space. 

We explore how the free parameters of the model affect 
the albedo and the climate of the planet. The planetary 
climate on Earth is most sensitive to variations in the liq- 
uid water cloud fraction. The albedo varies between 0.16 
and 0.58 for this parameter. The albedo varies between 
0.47 and 0.15, for precipitation efficiency between and 
100%. For number density between 70 and 300 cm -3 , the 
albedo varies between 0.25 and 0.47. The albedo varies 
between 0.26 and 0.3, for an ice cloud fraction between 
and 100%. The climate is least sensitive to changes in the 
relative humidity, unless it is below 20%. The albedo is 
between 0.3 and 0.27, for relative humidity between 20% 
and 100%. This effect is explained by the enhanced green 
house effect of high altitude clouds. 

Our microphysical cloud model provides a new and 
valid approach for modeling clouds and their effect on the 
climate for a wide set of rocky exoplanets within the hab- 
itable zone and will be used to study planetary science 
questions related to Earth. 
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